LPM effect as the origin of the jet fragmentation scaling in heavy ion collisions 
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We address a recent puzzling result from the LHC: the jet fragmentation functions measured 
in PbPb and pp collisions appear very similar in spite of a large medium-induced energy loss (we 
will call this "jet fragmentation scaling", JFS). To model the real-time non-perturbative effects 
in the propagation of a high energy jet through the strongly coupled QCD matter, we adopt an 
effective dimensionally reduced description in terms of the (1 + 1) quasi-Abelian Schwinger theory. 
This theory is exactly soluble at any value of the coupling and shares with QCD the properties of 
dynamical generation of "mesons" with a finite mass and the screening of "quark" charge that are 
crucial for describing the transition of the jet into hadrons. We find that this approach describes 
quite well the vacuum jet fragmentation in e + e~ annihilation at z > 0.2 at jet energies in the range 
of the LHC heavy ion measurements (z is the ratio of hadron and jet momenta). In QCD medium, 
we find that the JFS is reproduced if the mean free path A of the jet is short, A < 0.3 fm, which 
is in accord with the small shear viscosity inferred from the measurements of the collective flow. 
The JFS holds since at short mean free path the quantum interference (analogous to the Landau- 
Pomeranchuk-Migdal effect in QED) causes the produced mesons to have low momenta p ~ m, 
where m ~ 0.6 GeV is the typical meson mass. Meanwhile the induced jet energy loss at short mean 
free path is much larger than naively expected in string models. 
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Recently, the CMS Collaboration at the LHC reported 
the measurement of the fragmentation function of jets 
produced in Pb-Pb collisions at ^/s = 2.76 TeV per nu- 
cleon pair pQ. Surprisingly, in spite of the large energy 
loss signaled by the striking imbalance of jet transverse 
energies in the di-jet events [21 13] and by the suppres- 
sion of the large transverse momentum hadrons p|], it 
has been found that the jet fragmentation functions in 
Pb-Pb collisions are very similar to the ones measured in 
pp collisions. This is true even in the most central Pb-Pb 
collisions, and for the event classes where the imbalance 
of the jet energies is the largest pQ. It appears that the 
jet shape modification takes place only at small trans- 
verse momenta of the produced hadrons pt < 2 GeV, 
and at large angles outside of the jet cone pQ. This re- 
sult is surprising because the multiple scattering of the 
jet in the medium does not seem to produce extra high 
momentum hadrons. 

Let us briefly recall the space-time picture of jet evo- 
lution in QCD, see e.g. [5]. The produced high trans- 
verse momentum partons are in general far off mass shell 
and evolve emitting gluons and quark-antiquark pairs; 
this evolution towards smaller parton virtualities is gov- 
erned by the QCD renormalization group and described 
by DGLAP equations. At some scale Q\ ~ 1—3 GeV 2 the 
non-perturbative effects of confinement set in and trans- 
form the radiated partons into the observed hadrons. 
The uncertainty principle tells us that this "hadroniza- 
tion" occurs at longitudinal distance L/,. ~ zPj et /Q 2 , 
where z is the fraction of the jet transverse momentum 
Pj e t carried by the parton with virtuality Qq. Using for 
the sake of an estimate the values z = 0.2 , Pj et = 100 



GeV and Q Q = 2 GeV 2 , we get L h ~ 2 fm. In Pb- 
Pb collisions, this estimate suggests that the jet evolves 
down to the scales at which the dynamics becomes non- 
perturbative well within the produced medium. 

It is widely believed that the quark-gluon plasma at 
temperatures T < (2 — 3)T C produced at RHIC and 
LHC is non-perturbative and strongly coupled at scales 
of the order of 1 GeV. Combined with our estimate of 
the jet formation time, this suggests that to understand 
the CMS result one has to develop an approach to the jet 
interactions in the medium and its subsequent hadroniza- 
tion that is i) valid at strong coupling; ii) describes prop- 
erly the transformation of partons into the measured 
hadrons. Within the perturbation theory, the approach 
to the propagation of the jet in the medium has been 
developed in [SJ[7]; see OH] for reviews. It has been es- 
tablished that the Landau-Pomeranchuk-Migdal (LPM) 
effect [T0HI2] - the quantum interference of the radia- 
tion processes in the interactions with multiple scatter- 
ing centers in the medium - is as important in QCD as it 
was originally found in QED, although the non-Abelian 
effects modify the radiation pattern. While the LPM 
effect has been traditionally treated within the pertur- 
bation theory, it can be expected to affect the radiation 
amplitudes also at strong coupling. Indeed, the LPM 
effect may be viewed as a consequence of quantum me- 
chanics and the low-energy Low theorem that is based 
only on the symmetries of the theory (gauge invariance 
and the conservation of vector current) and is valid even 
when the perturbation theory does not apply. 

To proceed further we need an effective dynamical the- 
ory that can be treated at strong coupling and that can 
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describe the transformation of partons into the color- 
neutral hadrons. In spite of the lack of a complete theory 
of confinement much progress has been made in a qual- 
itative understanding of its possible key ingredients. In 
particular, it is likely that the "dual Meissner effect" pro- 
posed by 't Hooft [T3] and Mandelstam [TJ] describes cor- 
rectly the qualitative picture of the confinement of color 
electric flux. In this picture, the QCD vacuum contains a 
condensate of magnetic monopoles, and a quasi- Abelian 
Higgs phenomenon takes place. 

It has been suggested long time ago to use the (1 + 1) 
QED (the Schwinger model) [I5HT7] as an effective the- 
ory capable to model the transformation of partons into 
hadrons at high energies [TH] , see also [TH HHj ■ Indeed, 
the high energy justifies the dimensional reduction to 
(1 + 1) dimensions, and the Schwinger model with mass- 
less fermions shares many key properties with QCD, in- 
cluding: i) the Higgs phenomenon (local electric charge 
conservation is spontaneously broken); ii) the sponta- 
neous breaking of global chiral symmetry; iii) the screen- 
ing of color charge (similar to the scenario of confinement 
for QCD with light quarks proposed by Gribov [3T], see 
[22j for a review); iv) axial anomaly and the 0-vacuum. 
While the Schwinger model is Abelian, it may model the 
quasi- Abelian dynamics of QCD emerging due to the con- 
densation of magnetic monopoles in the vacuum that may 
be at the origin of confinement [T3J E] ■ The main ad- 
vantage of the Schwinger model with massless fermions 
(quarks) is that it is exactly soluble and can be used to in- 
vestigate e.g. the role of LPM effect at strong coupling. 
Of course, for our purposes we have to generalize it to 
allow for the rotation of quarks in color space as they 
traverse the quark-gluon plasma (recently the effects of 
the color flow were addressed in [23]). 

The Lagrangian of QED in 1 + 1 dimensions is given 

by 



where 



C = 



1 



(1) 



where is the U(l) gauge field, g is the coupling con- 
stant and if) is a Dirac spinor; x M = (t, z). The electro- 
magnetic vector current is given by j^(x) = ip{x)^ ll il}{x). 
It is well known that the theory ([!]) admits a bosonic rep- 
resentation in terms of a free massive scalar field theory 
with a scalar field <f> of mass m [THl HZ] : 
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In 



to describe e + e annihilation, an external current 



jg Xt of the produced quark and antiquark was added to 
the theory; the equation of motion becomes 



(□ + m 2 )(j)(x) = -m 2 4>ext(x) 



(3) 



(4) 



According to an external charge density corresponds 
to the scalar field 



4>ext{x) = vW dz'j° xt (t,z') 



(5) 



If we consider as an external source the pair of quark and 
antiquark moving along the light cone, the charge density 
is given by 

f ext (x) = S(z - t)6{z) - S(z + t)6(-z) (6) 

Given ([5]), the solution to ([3]), with boundary conditions 
<f>(t < 0, z) — and t 2 > z 2 , was shown to be [18] 

4>{x) =9(t + z)6(t-z)- A R (m 2 ,x 2 ) (7) 

where A R (m 2 ,x 2 ) is the retarded propagator. It is easy 
to check that the retarded propagator has the form 



A R (m 2 ,x 2 ) = J (m\x\)6{t + z)6{t - z) 



(8) 



where Jo is the Bessel function of the first kind and |x| = 
y/t 2 — z 2 . The solution to the equation of motion now 
reads 



4>(x) = 9(t + z)6{t - z)(l - J (m\x\)) 



(9) 



The (anti)kinks of the scalar field describe the production 
of quark-antiquark pairs in the vacuum induced by the 
external source. The bosonization allows us to describe 
this process in terms of the observed meson field <p. Let 
us evaluate the momentum distribution of the produced 
mesons. Eq. (|3| describes an interacting theory of the 
scalar field <f) with a classical source f{x). By using the 
Fourier representation of the free field and retarded prop- 
agator the momentum distribution of the created quanta 
is given by 
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where f(p) = J d 2 xf(x)e w ' x and E p = \/p 2 + m 2 . Let 
us now consider the quark and antiquark moving with a 
velocity v < 1 creating the charge density 



j°(x) = 5(z - vt)8(z) - S(z + vt)e(-z) 
Using ([5]) and ( [To| , we get 

r- —2vm 2 



f(p) 



E 2 — v 2 p 2 



from where 



2 4 
v m 



dN - 2 
dp ~ 71 E p (E 2 - v 2 p 2 ) 2 '' 
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note that there is no factor of N c here since we consider 
hadron yield per produced jet (with a fixed color ori- 
entation). For v = 1, we get dN/dp oc 1/E p familiar 
from the usual bremsstrahlung spectrum in (3 + 1) di- 
mensions. Let us define the usual fragmentation variable 
z = p h /pp et = p/pp et as the fraction of jet's momentum 
p 76 * carried by the hadron of momentum p h = p. As an 
application of the model, we use (p~3[) to evaluate dN/dz 



and fit it to the data on e + e — > hadrons at yfs = 201.7 
GeV [21]. The result is shown in Fig. [TJ the fit parame- 
ters are the scalar meson mass m and the matching scale 
Qo at which the DGLAP evolution has to be matched 
onto our model. From the fit we find m = 0.6 GeV that 
is consistent with the PDG value for the a meson, and 

Qo = 2 GeV (the velocity v = P jet / J 'Pj et + Q§). At 

small z < 0.1, our result is below the data points signal- 
ing the need for perturbative QCD evolution; however 
the spectrum at z > 0.15 is reproduced reasonably well. 



dN/dz 




FIG. 1: The spectrum of charged hadrons in e + e annihila- 
tion at y/s = 201 GeV; solid line is obtained from ( 13 1. 



Let us now extend this formalism to the case of a jet 
propagating through the quark-gluon matter. The quark 
will exchange color with the matter, rotating in color 
space and creating in the medium the static color sources 
located at coordinates Fig. [2J The different 

sectors in Fig. [2jare bounded by the sources with differ- 
ent orientations in color space; each of them is consid- 
ered as quasi- Abelian, and at large N c the production of 
mesons in each sector is independent. 

We see that there are only three different types of sec- 
tors in Fig. [2] - the one bounded by the quark escap- 
ing from the medium with no interactions (on the left, 
/1), the one bounded by the quark that underwent color 
rotation(s) in the medium (on the right, f^), and the 
one bounded on the sides by color static sources in the 
medium and a propagating quark from below (in the mid- 
dle, / 2 ). By using the methods described above, we get 
for the corresponding sources (we denote with t\ and ti 
the time when the first and second scatterings in medium 
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FIG. 2: The color flow in the jet interactions inside the quark- 
gluon medium. 
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Summing over the color orientations of different sectors 
(note that this does not bring in extra powers of N c in 
the 't Hooft limit of N c — > 00, g 2 N c fixed), we get the 
hadron spectrum 
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where we have omitted the interference between different 
sectors that is suppressed at large N c . We can write 
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In order to compare our result to the CMS data [T], we 
use the variable £ = In (1/z). Our result for the ratio of 
in-medium and vacuum fragmentation functions is shown 
in Fig. [3] We put TV equally spaced scatterings be- 
tween t\ and t% (the distance between them corresponds 
to the mean free path). One can see that the observed 
jet fragmentation scaling (JFS) is well reproduced if the 
mean free path \ m fp of the quark in the medium is short, 
A m /p < mT 1 ~ 0.3 fm. It is important to check whether 
the JFS in our computation stems from the absence of 
the energy loss - this would contradict the experimental 
observations [IHI]. The energy loss of the jet in medium 
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is large and consistent with experimental observations - 
this happens because the produced hadrons are abundant 
but possess small transverse momenta. Of course, our 
treatment has been to large extent model-dependent and 
can be both questioned and improved. However we feel 
that the origin of our result - the LPM effect at strong 
coupling and short mean free path - may be generic for 
a broad class of models. 
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FIG. 3: Ratio of in-medium to vacuum fragmentation func- 
tions. The length of the medium is fixed at 4 fm, the jet en- 
ergy is Ej e t = 100 GeV. Solid line: the first scattering occurs 
at t\ = 1 fm (assumed thermalization time), and subsequent 
scatterings occur with time spacing of At = 1/m = 0.3 fm. 
Dashed line: double scattering with ti = 2 fm and ti = 4 fm 
(At — 2 fm). Dot-dashed line: four scatterings with At = 1 
fm, ti = 1 fm. Data points are from pTJ. Open (filled) circles 
are for the leading (subleading) jet. 



is given by 



SE = 



dE h E h 



/ dN med dN vac 

V dE h dE^ 



[1 

[2 
[3 

K 

[5 

(17) [6 



We can use (17) to calculate the energy loss SE as a 
function of jet energy Ej et . We plot this in Fig. |2J note 
that our treatment is valid only when SE <C Ej et ; for 
short mean free path X m f p < 0.3 fm this means Ej et > 
100 GeV. The energy loss at X m f p < 1 fm is consistent 
with the values extracted from the data [T] , see [55] . 
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FIG. 4: Energy loss as a function of jet energy. The lines 
correspond to the parameters in the caption of Fig. [3] 



To summarize, we have evaluated the jet fragmentation 
function in QCD matter by using a non-perturbative ap- 
proach based on the exactly soluble quasi- Abelian model. 
We find that the observed jet fragmentation scaling is 
reproduced for any value of the mean free path that is 
shorter than 0.3 fm. Meanwhile, the induced energy loss 
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